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Possible Observations on a Signal ..... 


Patterns Used in Dual Parameter Model 


I. INTRODUCTION 


A. THE PROBLEM 

In several fields one must identify a true cause from a group of candidates. This 
identification must often be made with error prone information from a limited number 
of different observers or sensors. For example, in electronic warfare there may be 
several sensors which measure a signal from a radar, but exactly which radar type is 1t? 
In aviation, a pilot reports a rough running engine to aircraft maintenance. The 
mechanics use several testers to examine an engine’s performance and measure certain 
parameters. Because of poor calibration and operator error, very different readings 
result between testers. These errors may lead to “fixing what wasn’t broken and not 
fixing what was”. 

Let us consider our first fictious example, one occurring in electronic warfare. 
We are trying to identify a radar by the signal that was received. We can picture a 
problem in which we have J different radar types that may be in the area. Each of 
these radars 1s characterized by several parameters which can be measured and used as 
identification signatures. 

A radar has been detected. The parameters of that radar have been measured by 
I independent sensors. However, because of calibration errors, time since the last 
calibration, rough handling, operator error, atmospheric conditions, etc., the 
measurements do not seem consistent. We know, from historical data, that they tend 
to have a fair percentage of gross errors, or outliers. We do not know the specific 
error distribution that best represents an individual sensor at a given time. What 1s 
needed is a robust method of identifying the true signal, 1.e. one that can tolerate a few 
gross errors without being misled. 

A radar signal, received by observer 1 on parameter k, could be represented by 
the true signal 0, plus the error contributed by observer 1 namely Oy. The equation 1.1 


1 
would represent the value recorded by the cD observer on the kth parameter: 


Figure 1.1 shows a few possible outcomes from a problem with 3 observers and several 
radars. These samples represent the range of difficulty that might be encountered, 
from trivial identification to extremely difficult cases. In these examples radar #3 was 


actually active and the others were silent, but this fact was unknown to the observers. 


B. APPROACH 
1. Single Parameter Case 

If we assume that there is only one parameter available to identify the radars, 
then there are several Bayesian and Non-Bayesian algorithms that might be employed. 
We will compare a few of each for accuracy and robustness. We will consider various 
error distribution assumptions and look at the probability of each algorithm making a 
correct identification. The higher the probability of correct identification, or hit 
probability, the greater is the algorithm’s accuracy and usefulness. 

First, let us examine the single parameter case. Of course, our goal is to use 
analytical methods when possible. But, for the most part, simulation will be necessary 
to evaluate the following methods: 

Non-Bayesian Algorithms 
e average of the observations 
¢ weighted mean 
e standard median 
¢ voting algorithm 
e biweights 
¢ maximum likelihood method based on mixed normal distribution 
¢ maximum likelihood method based on Student - t distribution 
All of the above provide summaries of the ] determinations; the particular summary is 
then compared to each possible parameter (assumed known), and the nearest neighbor 
is identified. 
Bayesian Algorithms 
¢ based on normal distribution assumption 
¢ based on mixed normal distribution 
e based on Student-t distribution 
Each methods Bayesian posterior is computed for each parameter value, and the 


maximum posterior value is used to identify the active individual. 


1] 


we mea eam ee ee 


——— ey a me NS TT 


Va) i> 
= J a. 
<7 TJ 
»< x< 
< 
” x om C1 « 
< 
ra) ~~ om < N 
- Ee 4 —_— = 


ee 


poeb dar, Jd PRPIINNUDAOD LV AADNGOaMd|aY 


Possible Observations on a Signal. 


Preuiewlel 


If the assumed a priori probabilities are correct, then one would expect the 
Bayesian methods to supply better results. But, as the a priori- probabilities become 
less and less accurate these methods would be expected to be misled and become less 
effective. | 

In order to test the properties of the methods under consideration the error 
term will be generated using different distribution families. These error distributions 
will include: 

* normal 

e mixed normal with epsilons (mixing probabilities) equal to 0.1 and 0.25 

e Student-t with | degrees of freedom (1.e., a Cauchy) 
These different distributions have the affect of varying the weight and shape of the tail; 
the heavier the tails, the more outliers there are. The distributions and parameters 
used will be discussed more throughly in Chapter 2. 

Through simulations we want to examine the effects of the different error 
distributions on the identification methods. Also, we need to take a look at the 
sensitivity of the each technique to changes in the following: 

e scale (standard deviation for the normal model) 


e epsilon for the mixed normal based error model ( the probability of an outlier, 
roughly speaking) 


® prior probabilities for the Bayesian methods 
¢ number of observers (2 to 5) 
2. Bivariate Parameter Case 
A more realistic case to be considered is one in which there are two 
parameters, each of which can be measured and used to identify the radar. The 
measurement of these parameters will again be subjected to observation errors, but 
these errors may also be correlated. That is, the platform making the measurements 
may induce some error that 1s common to both parameters and some error that is not. 
We will use the standard correlation coefficient (p; ) to represent this interaction. 
We will examine the following four non-Bayesian algorithms for their 
robustness to outhers, changes in correlation, and scale interactions. 
e median 
e bivariate normal maximum likelihood estimator 
e bivariate mixed normal maximum likelihood estimator 


e bivariate Student t maximum likelihood estimator 
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Each of the above estimators will be evaluated using a mixed normal error distribution 
and a bivariate Student t, as in the one parameter case. These, as in the one parameter 


case, simulate a variety of conditions that we want to exafnine. 
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If. ERROR MODELS AND DATA GENERATION 


A. SINGLE PARAMETER CASE 
1. Model 

A radar emits a one parameter signal that is received and measured by I 
sensors. The | measurements of that signal are sent to a central processing point in 
order to identify the radar. Each measurement is an unbiased estimate of the true 
signal, which takes on one of J known values([), 5, ..., Hy). 

In order to simplify the analysis we have postponed treating the more realistic 
case of randomly spaced true signals and have concentrated on equally spaced Hi. 
Then, without loss of generality, we can assume that the radar signals are all | unit 
apart. Note that the computer code in Appendix B and C will handle the more general 
cases. For the Non-Bayesian algorithms the assumption that the signals are | unit 
apart creates a simple acceptance window or nearest neighbor decision rule. If the 
estimator of the signal (@) is within 0.5 of the true value, then a correct identification 
will be made. 

Though in actuality each sensor would have a distribution which models it the 
best, we have assumed that all the sensors have the same error distribution. This 
assumption allows for more direct comparisons between methods. This is not a 
constraint that 1s imposed by the computer code. The code in Appendix B permits 
different distributions for each observer. 

Another simplifying assumption was that the radar putting out the signals has 
no inherent variation. This is obviously not true but was reserved for later study. 

2. Data Generation 

Three families of distributions were used to simulate the error term in the 
measurements of a sensor. We will discuss each distribution in the following 
paragraphs. We will examine why these distribution were chosen and why the 
particular scale parameters were used. 

a. Normally Distributed Error 

The first distribution considered is the normal, with variances ranging from 
0.09 to 0.49. The normal distribution provides a strong base case for which to 


compare the estimators. It is very well understood, and in some cases the algorithm hit 
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probaoilities can be calculated or approximated without simulation, thus providing a 
means to verify our simulation procedures. 

The spread of variances from 0.09 to U.49, when combined with our 
previous assumption that the radur signals are all equally spaced at | unit apart Cowen 
a wide range of difficulties. For example. a variance of 0.25 puts 68:4 percent of the 
single observauons within our mentioned plus or munus 0.5 idenufication window. We 
would expect fairly accurate results. 

The computer simulation used the [MSL routine GGNMIL to generate the 
normal random numbers. The random numbers are all generated in subroutine 
DRAW. 

b. Alixed Normal Distributed Error 

The mixed normal combines two normal distributions. The first hus the 
same scale (standard deviation) as the previous normal; the second has a scale that is 
ten times as large. The second normal is chosen with probability € and the first with 
probability 1-¢;. The effect is that we induce a certain percentage of wide outliers. 
This ts another way of simulating heavy tails. Equanon 2.Iefsem Mood) Grambiiiyamd 


Boes [Ref. l:p. 122] was used to compute the mixed normal error model. 
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where: eae contamination factor 


scale of the first term 


aS 
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scale of the second term 


a 
~~ 
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Here the first scale parameter covered the same range as in the normal, but 
the second scale was 10 times the first. If the error contribution came from the second 
term of the mixed normal, then it 1s easy to see that a selection from 4 nOrinaleweeee 
standard deviation of 5 could result in a high probability of observations exceeding the 
range of our perceived signal values. Further, €, ranged from 0.1 to 0.25. This gives a 


wide range for the number of outliers one might expect: up to 25 percent. 
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Protea tne simmlatrom used the 1NISL routine GGNAMIL to generate the 
norma! random variables. GGUBIS was used to generate uniform random variable to 
pick whether the large scale or sinall scale parameter woulube used to compute the 
observations. 

c. Student-t Distributed Error 

Whe Siudemm® in-equation 2.2, was incitded to incorporate heavier tals 

than the normal. In this context the degrees of freedom were used as a method of 


adjusting the weight in the tails and 6, the spread. 


c(d.) (2.2) 
€ m = eee ae 
a ae (d+ 1)/2 

x.-8 l ] 

j ee 

1 - a oe 
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Where d. = degrees of freedom 
6, = scale parameter 
Values for the scale parameter were the same as 1n the normai case. It 1s 
. ae ° . . ° 
important to remember that the (scale parameter)~ 1s not the variance. The variance 1s 


given by: 
VAR(x:) = {d(d-2}6° for d.>2 (23) 


Note; this equation does not apply for degrees of freedom 2 or less. The Cauchy, 
which was used as one of our error sources, is a t with degrees of freedom equal to | 
and has an undefined “infinite” variance. See Barr and Zehna A Cauchy has 
notoriously heavy tails and should pose a difficult test for the algorithms. 

The computer model-in subroutine DRAW of Appendix B uses the IMSL 


routine GGAMR to generate the pseudo random numbers. 


B. BIVARIATE PARAMETER CASE 
1. The Model 
A radar again emits a signal. This signal contains two parameters which are 


; : 2 ») ' 
measured, Each estimate of a parameter has some variance 6~,, or 6~,,, associated 
“| 


I 
with it; wnere the index 1 refers to the jth observer. There is also some covuriance that 


relates the interaction between the two estimates. We will use the classical correlation 


coefficient to represent this interaction. 
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There are two descriptions of the placement of the pairs that describe a radar. 
These are shown in figure 2.1 
We will maintain only one error distribution for all observers. not because we 
bel:eve this to be the case, but because it makes comparison easier. 
2. Data Generation | | 
a. Bivariate Normal Error and Mixed Normal Error 
The bivariate nuxed normal distribution, given in equation 2.4, will be used 


as one representauon of the errors. 
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Tne equation 2.5 from Law and Kelton [Ref. 3:p. 269] was used to sSimiiiaieuae 


observations. 
x,” N(@ ,o a 
j a x de 
i (2.5) 

ey 

l 2 2 
S ee -6 10,0 1- 
yi Js, + Ca eat N[ : (1-p,)] 
x, i 


where: ex, and Sy. are the scales associated with the first or second terms 
The values for epsilon range from 0 ( bivariate normal ) to 0.25. The py; (correlation 
coefficient for the first term) and py, (correlation coefficient for the second term) will 
France 11 orn =o vers. 
b. Bivariate Student T Error 
The Bivariate Student T, as in the single parameter case, poses a difficult 


test for the algorithms. 
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PATTERN 1 





PATTERN 2 





Figune 2.1 


Patterns Used in Dual Parameter Model. 
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Equation 2.7 from Kotz and Johnson [Ref. 4:p. 129] was used to simulate the 


observations. 
Py 4 ry =~ 
x= Z.(Yid A) O + Hi oe 


h observation from a bivariate normal 


Where: Z; 1s the it 
Y 1s distributed chi-squared with d; degrees of freedom 


The values for p range from 0.5 to -0.5. 


Il. EXANIENATION OF METHODS 


A. SINGLE PARAMETER CASE 
1. Non-Bayesian 
a. The Nlean 
The mean-nearest neighbor algorithm was incorporated into the studv to 
examine how it compares to what we believe to be more robust estimators. Also, we 
can calculate the probability of a correct identification given normally distributed and 
mixed normally distributed error analvtically. The mean considered here is a simple 
average of the observations and is done in the subroutine labeled MEAN. 
Mf We assume the error is [rom a normal distribution, then we can calcuiate 
the probability of a correct answer for a specified number of observers using equation 
Bel. 


1(7- oC 


By +$/2 _ 7 
Pal els oe 3 (3.1) 
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where ¢= Bs aa 


Oo = standard deviation of the mean 


Then, using standard notation, the probability of a correct identification can be 


calculated by equation 3.2. 
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We can also solve for the EeGHE) Ea \ienmenecmeniotuis assumed to be 


mixed normal by using equation 3.3 (assuming €,= €, 1= 1,2,...,I) 
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The results for an interior point where 6, is 0.5, a 6» is 5 and when I is ranged from 2 
to 20 is contained in Table 1. Note that as the number of observers goes from 2 to 6 
that the probability of hit goes down! This trend reverses as the number of observers 
gets larger than 6. This is a result of the large outliers effect on the mean, as modelled 
by the € contamination. In reflects the unfortunate non-robustness of the 
mean-nearest neighbor algorithm. The probability of getting a draw from the side of 
the mixed normal with a large scale increases with the number of observers. The mean 
can not overcome these big outliers with a limited number of observers, but eventually 
the central limit theorem effect takes over. This example strongly suggests that 
ordinary linear (mean-based) rules will be ineffective in the presence of outliers. The 


algorithm to compute the mean is contained in subroutine MEAN of Appendix B. 


TABLE 1 
Pry FOR MEAN WITH MIXED NORMAL ERROR 
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b. The Weighted Mean 
The weighted mean is identical to the average when the variances of all the 
observers are thought to be equal. But if the variances are believed to be different, 
then the observation from the watchers with the lower variance should carry more 
weight than those with a higher variance. Equation 3.4 is the appropriate weighted 


mean. 
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c. The \Iedian 
The standard median was also introduced, as a simple robust estimator. 
The median has the advantage of requiring no prior knowledge of the sensors or the 
radars. Furthermore, we can approximate the probability of a correct guess. If we 
assume normally distributed error, then from Ilodges and Lehman [Ref. 3:pp. 926-927] 


and Chu [Ref. 6:pp. 112-116] we get equation 3.5. 


2 
oC 
Ze J 1 aS 
Median NH. Tit aaa ( ) 
1 I 
4 ; 
—- 1 if I is odd 
TT 
where a = 
6 
| 7h if I is even 


The median is found by the computer function RMEDIAN, but 1s called by subroutine 
MEAN, ; 
d. Voting Algorithm 

The voting algorithm is a democratic process. For the present uniform 
parameter distribution this process picks the radar that has the most observations 
within 0.5 units of it. If there is no one radar that has more that the others than the 
median is chosen. 

It is easy to see that this results in the same choice of a radar as the 
median for the two and three observer cases. The only advantage is when there are 
four or more observers. Then the ability of the voting algorithm is not to be mislead 
by two outlier observations that happen to be on the same side of the true value. On 
the other hand, any two observations(even if outhers) that are close to a b will, more 


than likely result in that radar being identified as the one. 


e. Biweights 


The biweight procedure 1s a robust iterative process that uses the calculated 
A : ° . . « « « 
Q‘(t) to begin iteration t+ 1. Our first guess is at the median of the observations. The 
updating process proceeds as in equation 3.6. 
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The constant c was set equal to 6, as recommended in Mosteller and Tukev (Ref. 7:p. 
353] This algorithm is contained in subroutine BIWGHT of the simulation in Appendix 
B. 

ff. Maximum Likelihood Estimator Based on Mixed Normal 


[f we anticipate that all the observations have a mixed normal error 


distribution as described by equation 2.2;7~then the maximum likelthood method 


described in Mood, Graybill and Boes [Ref. I:pp. 276-286] proceeds by equation 3.7. 
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Next We take the natural logarithm of cquation 3.7 and rearrange terms resulting in 
equation 3.8. 
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Differentiating both sides we get equation 3.9. 
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We apply the iterative reweighting method simular to that used in the biweights to get 


equation 3.10. 
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‘As in the biweights, we again use the median of the observations as the starting point 
of our weighting. 

iivemae oc mermmodmisecanieca out in subroutine EPSMLE. The process 
should converge and stop. For computer purposes, convergence is considered to be 
when 09(t)=@(t+1)+-.0001. If the process did not converge by the 100g iteration, 
then 1t would stop and take the median. Also, a flag would be sent up to warn the 
operator that this run 1s not going to be as completely a result of the expected method. 

g. Maximum Likelihood Estimator Based on Student t 
Tie wert edsed on the Student t assumes that each observation has a 


Student t error distribution. To solve the MLE we do the following: 
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Now tuking the natural log we arrive at cquation 3.12, after discarding inpeleaee 
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Unfortunately, this could have several solutions. We can usually get an appropriate 
solution using iterative reweighting. Again, the startuung point will be the median of the 
observations. Note that if 0(r) is replaced by the median of the observations, there will 


always be a unique solution. 
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This procedure is carried out in subroutine TMLE. This method, as in the one above, 
stops when @'(t)=0’(t+1)+-.0001. If, after 1000 iterations, it has not converged, then 


1t selected the median and set a flag as before. 
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h. Identification 
In ali the Non-Bavesian methods the resulting value was an estimiute of the 
true value. To state which radar type we actually heard a comparison was made. Tius 
involved picking the radar t\pe wiose true Value was closest to the estimated value. If 
fe estiimatca Valuéeexcecded the largest H; or was less than the smallest Me. then the 
identified j would be the one on the appropriate tail. This is the nearest neighbor 
procedure. 
2. Bayesian methods 
a. Where x; is Normally Distributed 
ere tia Oceacutttonal iformation that we want to take advantage of 
which was not incorporated in the preceding methods. That is, there might be assumed 
Knowledge of the probability of a particular rudar type being in the range of-detection. 
Pomecwmpien.c May know that half of the radars in the area are of one particula 
type. Then, all else being equal. we would expect to gain detection of that particular 
type more often than on the others. We can use this information as the a priori 
probability for a Bavesian approach. With normally distributed observations the Baves 


posterior 1s: 
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We can simplify to get equation 3.17. 
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The bavesian normal method 1s contained in subroutine NORMBY in Appendix B. 
b. Where x; has a Mixed Normal Distribution 
If we contaminated the above error by adding, with some probability 
epsilon, wide outliers from a mixed normal distribution: then we can again solve for the 


posterior shown in equation 3.158. 








I a z 
P | 
p..(6 ) = IKp oe (1-¢ , ) | exp | 3/3 4 | (3 S) 
i=] vine Oey il 
Ne eXp e 
om op 1zZ 








* p.m 
ies 24 
OY). = gay 19) 
p, (8) ; (SPE 
m. 
> Pym 
j=1 
where 
I 
ifs 
m. = ex - LN(S % a 
j i > meee Cie) sa : ee 
i=1 
1 2i 1 1 
Ce. see .) — 
nee eG | 
+ LNieé. + (1-<€ ae, LE 
1 * Gece i1 © 42 


il 


The Bavesian mixed normal method is contained in subroutine EPSBAY in Appendix 
B. | 
c. Where x; has a Student t Distribution 
We can follow the above arguments with a t-distributed error as well. 
Where equation 2.2 1s the underlying model of the error. We substitute @imegeage 


equation for the posterior and get equation 3.20. 
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This routine is performed in subroutine TBAYES in Appendix A. 


B. TWO PARAMETER CASE 
[. The Median Based Algorithm 

Many of the algorithms used in the single parameter case can be generalized 
to more dimensions. We have limited ourselves to only four algorithms to be examined 
in the tWo parameter case. All four are non-bavesian and will be compared by their 
probability of correctly identifving a 2 parameter radar. 

The first algorithm is a method based on the median. We proceed by finding 
the median of the I observations on the first parameter and the median of the 
observations on the second parameter. The euclidean distance from this pair to all the 
possible pairs is determined. Then we take the radar that is closest to the medians (the 
nearest neighbor) to be the radar that is sending. 

This method is carried out in subroutine RMED. 

2. Bivariate Normal Maximum Likelihood Estimator 

If we assume the errors to have a bivariate normal distribution with known 

parameters, then we can calculate the likelihood of cach parameter pair. This can be 


maximized to determine an estimate of the value in the usual way as shown in equation 
oD. 
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This procedure is accomplished in subroutine N VIVEZ. 
3. Bivariate Mixed Normal Maximum Likelihood Estimator 
One would not expect the above method to be very robust. We can improve 
on this bv expanding the algorithms of the single parameter case to the two parameter 
case. Equation 3.25 shows an MLE method based on the bivariate nuxed normal 


distribution. 
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This algorithm 1s contained in subroutine EPMLE2 in Appendix C. 
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4. Bivariate Student t Nlaximum Likelihood Estimator 
Another algorithm can be created based on the bivariate t MLE. Using the 
bivariate Student t in equation 2.6 we can follow the basic mle method to construct the 


new estimator in equation 3.24. 
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This procedure is carried out in subroutine TMLE2 in Appendix C. 
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IV. CONCLUSION AND RECOMMENDATIONS FOR FUTURE 
RESEARCH 


1. Observations 

We examined several identification methods for the single and double 
parameter cases. Each algorithm was tested under various conditions of true and 
assumed (value used in the algorithms) scale, priors in the single parameter case, and 
pattern and correlation, p, in the two parameter case. 

In the single parameter case we examined three Bayesian and seven 
non-Bayesian methods. The Bayesian methods included one based on a normal, one 
based on a mixed normal (with two different values for epsilon) and one based onat 
(with three different values for degrees of freedom). The non-Bavesian methods 
included a simple average of the observations, a variance weighted mean, median of the 
observations, voting algorithm, biweight and MLE’s based on the mixed normal and 
Student’s t. As expected the Bayesian algorithms did slightly better than their 
non-Bayesian counterparts if the priors assumed by the Bayesian were correct. But, if 
the priors were not correct then the Bayesian methods were significantly mislead. 
Therefore, the Bayesian algorithms were not reconimmended unless reliable priors are 
available or vague priors are used. The algorithms all displayed a tolerance to 
Inaccurate scale parameters with only a 1 percent decline in the probability of 
identification. The mixed normal and Student’s t based method (which were designed 
to handle outliers) did as well as the average and median on the normal based 
Observations and significantly out performed them when the number of outliers was 
increased with the mixed normal and t error distributions. Therefore, for the single 
parameter case either the mixed normal or Student’s t based MLE methods are the 
FeCOMMmended mROcedunes. 

In the double parameter case we concentrated on the non-Bayesian methods. 
These methods included an algorithm based on the median and three MLE’s; one 
assuming bivariate normal error, one assuming bivariate mixed normial error and the 
lasted assuming bivariate Student’s t error. The probability of a correct identification 
was very dependent on the true correlation coefficient and pattern of the observations, 
in addition to the affects of the other elements in the single parameter case. As an 


example, a linear pattern in combination with a negative correlation, p, resulted in a 9 


Ga 
tJ 


percent increase over the box pattern. But, a positive p in those conditions came up 
with 4 percent loss for the diagonal pattern compared to the box pattern. The 
algorithms adjusted well to inaccurate information, except for outhers. The normal mle 
and to a lesser extend the median could not cope with the outliers, when compared to 
the t or mixed normal based mle’s. As in the single parameter case, the t or mixed 
normal based mle is strongly recommended as the algorithm of choice. 

2. Areas for Further Research 

There may be more information available to the user from these algorithms 
than just which radar is thought to be radiating. For the Bayesian there is a 
probability associated with each possible radar. We can also determine a posterior for 
the MLE methods. If we divide the value that results from solving the algorithm with 
the known parameters values in it, by the sum of all these values for all the radars, 
this gives a posterior equivalent to the Bayes if all a priori probabilities were equal, i.e., 
for a vague prior. These two methods can be used to indicate how confident we should 
be in a particular choice. 

Another possiblity for the non-Bayesian estimates is the absolute deviation 
from the estimate of the parameter to each of the radars. For example, if the estimator 
was on the border line we should be less sure about identification than if it is directly 
on. Another method might be to remove the observation that is furthest away from 
the non Bayesian estimate. Then recompute the estimate based on one less 
observation. Do the results change significantly? 

JihieSeWproeedures can also™@ibe peneralized to include up to p (p> 2) 
parameters, while also permitting parameters to be missing from from the results of 
some of the observers. For example; observer 1 might record observations on the first 
and third parameters, observer two on the second an third and observer 3 on all of 
them. This case is not currently studied, but a liklihood-based analysis seems possible. 

The affects of variance in the signal send by the radar also needs should be 


looked at. How does this affect the estimators? Are some more resistant than others? 
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APPENDIX A 
SIMULATION RESULTS 


I. SINGLE PARAMETER CASE 

The following tables show the results for a wide range of possible cases that were 
investigated. The tables represent the outcomes for the probability of a particular 
algorithm correctly identifying the radar that is sending. Each entry is the result of 10 
runs of 1000 trials each. The entries are expected to be within +- .01 of the true 
values. The tables show the results for the four error distributions: normal, mixed 
normal (epsilon=.1), mixed normal (epsilon=.25), and Cauchy. These are listed 
across the top pf the tables. Also, on the top are the number of observers used to both 
generate and analyze the data. Down the left side of each table are the algorithms and 
their associated degrees of freedom or epsilons. Those above the double lines are 
Bayesian and those below are non-Bavesian. More details and results from each table 
are presented in the following paragraphs. 

Tables 2, 3 and 4 show the results from the best case, where the parameters used 
to analyze the observations are exactly correct. The only time this 1s not true is whcre, 
in the t based estimators or the mixed normal based estimators, the degrees of freedom 
or epsilon have been varied. These changes are represented in the left hand column. 
The priors for these tables were all equally likely. An observation was produced by 
sclecting a mean, based on the priors, from those given; then generating the error from 
the appropriate distribution and finally adding the two. Each draw for a mean and 
error are independent of all the others. 

Under ideal circumstances, the Bayesian methods have a slight advantage over 
the non-Bayesian procedures for larger numbers of observers, and a very significant 
edge for the lower number of observers. This was expected because of the additional 
information available to the Baycsian estimators. What is of more interest is the 
sensitivity of the algorithms to inaccuracies. 

In tables 5 and 6 the sensitivity of the methods is compared. Here the error 1s 
generated with a variance or spread of 0.25 and the estimators assumed 0.09 or 0.49 
respectively. The priors are all equal and the mean of the observation (radar in action) 


is Selected at random from the 5 given. 
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We see that all the algorithms are fairly insensitive to changes in scale. But, of 
importance is the apparent lack of change in both the t and median algorithms. While 
the normal and mixed normal show a 2 to 3 point degradation. 

In tables 7 and 8 we examine the effects of varying the scale parameters. That is 
to say, the true scale for the first receiver was 0.4, the scale for the second was 0.5, the 
scale for the third was 0.6, the fourth was 0.7 and the fifth was 0.8. The results in 
Table 7 show the probability of a correct identification when the correct scales are 
known and used in the algorithms. Table 8 shows the results if the correct scales are 
not known, but the average of the scale parameters 1s used. The difference 1s small. 

Tables 9, 10, 11 display the results from simulation when the priors for the radars 
doing the announcing are inaccurate in various ways. For these radar 3 was always 
sending, but in table 9 the priors assumed by the Bayesian estimators were 0.15 for the 
first, 0.3 for the second, 0.1 for the third, 0.3 for the fourth, and 0.15 for the fifth. In 
Table 10 we were less accurate with the first having a prior probability of 0.3 , the 
Saeomame! O11, tie third of 0. the fourth of 0.15 and the fifth of 0.3. Finally, we 
introduced a set of priors probabilities that was not symmetric with the first as .05, the 
second as 0.1, the third as 0.15, the fourth as 0.3 and the fifth as 0.4. 

Changes in the priors created the largest degradation in the probabilities of 
making a correct identification. The Bayes procedures are extremely sensitive to 
inaccuracies in the priors. They lose the significant advantage they held in the two 
Observer cases. The priors tend to draw away the Bayesian methods when the true 
signal is close to another signal that has a larger prior. This represents a particularly 
significant distraction to the Bayesian methods. 

Tables 12 and 13 display the probabilities of correct identification if radar 1 1s 
always sending or if radar 3 1s always sending and the priors are assumed to be equallv 
likely. 
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Cn 


ALGORITHM 


BAYES NORM 


BAYES EP.1 


BAYES £.25 


BAYES T1DF 


BAYES T3DF 


BAYES T10D 


AVERAGE 


WEIGH MEAN . 


MEDIAN 


TABLE 2 


Py WITH SCALE PARAMETERS OF .3 AND PRIORS EQUAL 


ERROR DISTRIBUTION 
NORMAL EPSILON .1 EPSILON .25 STUDENT T 1DF 


5/2 573 5/4 5/5 5/72 573 574 5/5 5/72 573 574 5/5 5/72 5/3 5/4 5/5 


98 1 = 1 .88 .86 .85 .85 -72 .70 .68 .68 72 «735 .72 See 
=o 868 : 1 94 290 2 99 oe 66 .93 .97 .98 -§60 .88 .91 .9% 
Ie OL 1 1 oo .9o 2.99 1 J00 .9D «97 290 .80 .88 .91 .93 
96" 1 Z L 94 296 .99 1 65 .95 .9e .95 81 .90 .95 236 
mo . A oS See Fo -O9 .9e «MO oF 81 .90 .935 Je 
vor 1 : aL We 69 IO ID 80 .91 .9% .96 79 .B89 .91 98 
98 1 Z J 88 .87 .85. 2.86 of > «AD» OB 69 735.735 «135 «oe 
25 2 i 88 .87 .85 J65 73 27 ORG OMG9 ~73 373 75a 
98a 99 ot a #88: 697 «98 .99 «73 Ieee 73 .88 .89 .9¢ 
“985.99 ot 1, 8 .W7on98. .97 «MD » In sBbee96 -73 .88 .90 .94 


VOTING 


BIWEIGHT 


MLE T 1DF 


MLE T 3DF 


MLE BPS .1 


MLE EPS.25 


NOTE: 


OM PWD eH 


no6. 99 1 88 .96 .98 .99 73 (CB) .985..95 73 .86 .89 .98 


98.99 1 1 .88 .98.99 1  .73 .92 .94 .97 .74 .89 .91 .95 
98 1 od 1 BB . 97s 9SecDecce7Suu9]) 996). 76). Seuenee 
98 1 1 1  .88 .98 .99 1 .73 .92 .94.98 .73 .89 .90 .94 
98 1 1 1  #.88.98 .99 1 .73 .92 .94 .98 .73 .88 .90 .% 


PRIOR PROBABILITIES ARE ALL EQUAL 

SCALE PARAMETER WAS .3 FOR ALL ERROR DISTRIBUTIONS 

FOR MIXED NORMAL ERROR, 2nd SCALE WAS 3 

THE ASSUMED SCALE WAS EQUAL TO THE TRUE SCALE 

ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
ALL ENTRIES WERE ROUNDED OFF TO 2 DIGITS 
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TABLE 3 


Pi WITH SCALE PARAMETERS OF .5 ANO PRIORS EQUAL 


ERROR DISTRIBUTION 

NORMAL EPSILON .1 ERSELON...25 STUDENT T 10F 
ALGORITHM 5/72 5/73 5/74 5/75 5/2 5/73 5/74 5/5 5/2 5/73 5/4 5/75 5/2 573 574 5/75 
BAYES NORM .88 .93 .96 .98 077 «6/8 «178 OST .64 .60 .58 .57 -60 .59 .60 .60 
BAYES EP.1 .88 .93 .96 .98 82 .90 .9% .96 -7> aes «OS .92 som ./5 2/9 .83 
BAVES E.25 .87 .93 .96 .98 .82 .89 .94 .96 .75 «85. 766 «92 sO7ect1oe.79 83 
BAYES .T1IDF .86 .92 .95 .97 -6l wes .93 .95 75 .82 .87 .91 -68 .77 .82 .86 
BAYES TSOF .88 .93 .96 .98 sCeuee? 295 «95 74 .82 .86 .90 J64) «764.61 .35 


BAYES T10D .88 .93 .96 .98 .80 .88 .92 .9% 70 .78 .83 .87 -65 .74..79 .83 


AVERAGE -88 .93 .96 .98 off 0 fB 6179 377 -6% .61 .59 .58 .60 .60 .61 .60 
WEIGH MEAN .88 .93 .96 .98 7? .78 .79 .77. =.6% .61 .59 .58 .60 .60 .61 .60 
MEDIAN 88 .89 .95 .95 fies «9 .92 .64¢ .78 .81 .86 .60 .75 .78 .83 
VOTING 88 .89 .91 .95 7/7 .85 .89 .91 .6% .78 .80 .85 .60 .75 .78 .83 
BIWEIGHT 88 .90 .9%4 .96 .77 .84 .91 .93 -64 .76 .81 .85 .60 .73 .77 .8l1 
MLE T 1DF .88& .91 .95 .97 77 .88 .92 .95 .64 .81 .85 .90 -61 .76 .80 .85 
MLE T 3DF .88& .93 .96 .98 .77 .88 .92 .95 .65 .80 .8¢ .89 .61 .76 .80 .8% 
MEE EPS .1 .88 .93 .96 .9%8 7789 .935m. 96 .64 .82 .86 .91 .60 .74 .78 .82 
MLE EPS.25 .88 .93 .96 .98 77 .89 .93 .96 .69¢ .82 .86 .91 .60 .75 .78 .83 
NOTE: 

PRIOR PROBABILITIES ARE ALL EQUAL 

; SCALE PARAMETER WAS .5S FOR ALL ERROR DISTRIBUTIONS 

FOR MIXED NORMAL ERROR, 2nd SCALE PARAMETER WAS 5 

FOR T OISTRIBUTEO ERROR, DEGREES OF FREEOOM WERE 1 

THE ASSUMED SCALE EQUALEO THE TRUE SCALE 


ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
ALL ENTRIES WERE ROUNOEO OFF TO 2 OIGITS 


NOU PWD 
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TABLE 4 


Py WITH SCALE PARAMETERS OF .7 AND PRIORS EQUAL 


ERROR DISTRIBUTION 

NORMAL EPSILON .1 EPSILON .25 STUDENT T 1DF 
ALGORITHM 5/2 5/73 574 5/75 5/2 573°574 575 5/2 573 5/746 5/5 5/2 5/73..5/7@ 5/75 
BAYES NORM .75 .83 .87 .91 .66 .69 .70 .71 055 «52 352 .50 -52 .51..51 oe 
BAYES EP.1 .75 .82 .87 .91 71 .78 .84 .88& 69 .72 .77 wee 58 .63 .68 .73 
BAYES E.25 .75 .S25.57 0 -71 .78 .83 .88 64 .72 .78 .82 58 .635...69 OS 
BAYES TIDF .73 .80 .85 .88& -69 .76 .82 366 -63 .71 .76 .80 60 .66 .72 .77 
BAYES TSDF .75 .81 .87 .90 71 .77 .83 .87 64 .71 .76 .80 58 .66 .71 .76 


BAYES TLOD .75 .82 387 .& .69 .76 .81 .86 -60 .67 .73 .77 55 .63 .68 [a 


AVERAGE a> . osaxe87. . 9 -66 .69 .70 .71 56 .53 w52 .51 52 .52 Siege 
WEIGH MEAN .75 .83 .87 .91 -66 .69 .70 .71 56 .53 .52 .51 52 .52 .51,.52 
MEDIAN .75 .77 .84 .86 -66 .73 .79 .82 56 67 71 .74 52 .65 .68 .74 
VOTING 75 .77 .79 .85 -66 .73 .75 .8l1 -56 .67 .68 .72 -52 .65 .67 .73 
BIWEIGHT 75 .78 .84 .88 .66 .73 .81 .85 -56 .66 .72 .75 52 .63 .67 .72 
MLE T 1DF .75 .79 .85 .88& -67 .75°°6. .&5 -56 .69 .74 .79 53 .66 .71 Fae 
MLE T 3DF .75 .81 .87 .90 .67 .77 .82 .86 -56 .70 .74 .79 -53 .65 .70 .76 
MLE EPS .1 .75 .82 .87 .91 -66 .78 .83 .88 56 .71 .76 .82 52 .63 .68 .73 
MLE EPS.25...75°.82°.87" .90 -66 .78 .83 .88 -56 .72 .76 .82 52 .64 .68 .73 
NOTE: 

PRIOR PROBABILITIES ARE ALL EQUAL 

SCALE PARAMETER WAS .7 FOR ALL ERROR DISTRIBUTIONS 

FOR MIXED NORMAL ERROR, 2ND SCALE PARAMETER 49 

FOR T DISTRIBUTED ERROR, DEGREES OF FREEDOM WERE 1 

ASSUMED SCALE EQUALED THE TRUE SCALE 


ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
ALL ENTRIES WERE ROUNDED OFF TO 2 DIGITS 
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TABLE 5 


Pu WITH ASSUMED SCALE OF .3 AND TRUE SCALE OF .5 


ERROR DISTRIBUTION 
NORMAL EPSILON .1 EPSILON .25 STUDENT T 10F 
ALGORITHM 5/72 5/73 5/74 5/5 5/2 5/3 5/74 5/5 5/2 5/73 574 5/75 5/2 5/73 574 5/75 
BAYES NORM .88 .94 .96 .98 we’ .7/8 .78 217 woo «59 m5? 255 Poo. 57 9 


BAYES EP.1 .84¢ .90 .93 .95 .80 .87 .91 .94 -f> «61 .86 ~.90 65 .74 .79 .82 


BAYES T1DF .85 .91 .94 .96 -80 .88 .92 .95 To eOL .O7 o9L -67 .77 .82 .86 
BaYES@RSDF 967 92 395 [97 we2 88S 1693 095 7% .82 .87 .91 68 .77 .83 .86 


BAMES T10D 287 «93 «96 «9/7 aBZenS9 «9S «96 -74 .81 .87 .90 -67 .76 .81 .85 


AVERAGE -88 .94¢ .96 .98 wf! 70 (GS at? -64¢ .61 .59 .58 -60 .60 .60 .60 
WEIGH MEAN .88 .94¢ .96 .98 047 78 78 OT 64 .61 .59 .58 -60 .60 .60 .60 
MEDIAN woo 609 395.95 ary 6 ©OpmeIO 692 -64 .78 .82 .87 -60 .75 .78 .83 
VOTING -66 .S9 .92 .95 Wino > .S6".92 -6%¢ .78 .81 .86 -60 .75 .78 .82 
BIWEIGHT 88 .90 .9%4 .96 -/7 .85 .91 .9% .64 .76 .82 .86 60,.7c af? .S1 
MLE T 10F .88 .90 .9% .96 -78 .87 -..92 .94 .6¢ .80 .85 .89 -61 .76 .81 .85 
MLE T 30F .88 .92 .96 .97 -/8 .88 .93 .95 -64 .80 .85 .90 -61 .76 .81 .85 
MLE EPS .1 .88 .89 .93 .95 wi? guy 70. oS -64 .80 .83 .89 -60 .75 .78 .83 


BAYES E«25 .835 .89 292 295 279 @B6 .90 .93 7/2 .80 .85 .89 64 .73 .78 .81 
| 
| 

MEENErowc> soo .6G .92 .93 ,f?7 «66 6.69 . 92 -64¢ .79 .81 .88 -60 .74 .77 .82 | 

| 


NOTE: 

PRIOR PROBABILITIES ARE ALL EQUAL 

TRUE SCALE PARAMETER WAS .5, FOR ALL ERROR DISTRIBUTIONS 
FOR MIXED NORMAL ERROR, 2ND SCALE PARAMETER WAS 25 

FOR T DISTRIBUTED ERROR, DEGREES OF FREEDOM WERE 1 
ASSUMED SCALE WAS .3 FOR ALL ALGORITHMS 

FOR MIXED NORMAL ASSUMED 2ND SCALE WAS 3 

ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
ALL ENTRIES WERE ROUNDED OFF TO 2 DIGITS 
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ALGORITHM 


TABLE 6 


Pi WITH ASSUMED SCALE OF .7 AND TRUE SCALE OF .5 


ERROR OISTRIBUTION. 


NORMAL EPSILON .1 EPSILON .25 STUDENT J 10F 


5/2 5/73 5/74 5/5 5/2 5/3 5/4 5/5 5/72 5/73 5/4 5/5 # £572 5/73 5/4 5/75 


BAYES NORM .88 .93 .96 .98 77 .79 .78 .78 #3 .62 .59 .58 .60 


BAYES EP.1 .68 .93 .96 .98 -82 .89 .93 .% 73 .61 .65 769 . 66 


BAYES E.25 .88 .93 .96 .98 -82 #89 298 .96 43 .82 .87 .91 -66 


BAYES T1OF 


BAYES T3DF 


87 .92 .95 .97 <Gey.69).95 «95 -74 .82 .87 .91 68 


#68 .93 .96 @6 ie) .89 .92 .95 72 .B6Y .8B .89 66 


BAYES T100 .88 .93 .96 .98 -78 .86 .90 .93 65.75 29.65 63 


-61 .60 .60 


7% .77 .81 


~-f5 ./60e2 


.7/8 .81 .86 


./7 .80 .85 


AVERAGE 


88 .93 .96 .98 47.19 tee. 1G -63 .62 959 .56 .60 


WEIGH MEAN .88 .93 .96 .98 47 . 79 =—78 ./8 -65 .62 <59 .58 60 


MEDIAN 
VOTING 
BIWEIGHT 
“MLE T 1OF 


MLE T 


30F 


88 .89°%9G%95 4W77Meo Fol o2 Woz emteate. Wen 
.88 .89 .91 .95 .77 .86 .89 .92 .63 .78 .81 .85 .60 
.88 .90 .94 .96 .77 .85 .91 .93 .63 .77 .82 .86 .60 
.88 .92 .95 .97 .77 .89 .92 195 "Wes veil es 90 Wer 


68 .95 .96 .98 ./7 .88 «92 .95 .64¢ .80 .84¢ .88 61 


MLE EPS .1 .88 .93 .96 .98 «ft? .89 «92575 -63 .80 .84 .89 60 


NEEVEPS.25 .66™.95 <96 398 -77 .89°.93 .96 .63 .81 .86 .90 60 


NOTE: 


On AML WP 


PRIOR PROBABILITIES ARE ALL EQUAL 

TRUE SCALE PARAMETER WAS .5, FOR ALL ERROR OISTRIBUTIONS 
FOR MIXED NORMAL ERROR, 2ND SCALE PARAMETER WAS 25 

FOR T DISTRIBUTED ERROR, DEGREES OF FREEDOM WERE 1 
ASSUMED SCALE WAS .7 FOR THE ALGORITHMS 

FOR MIXED NORMAL MLE 2ND SCALE WAS 7 

ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
ALL ENTRIES WERE ROUNDED OFF 170 2 OIGITS 


40 


-76 .77 .84 
76 ./77 .83 
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Pu 
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ALGORITHM 5/72 5/73 574 5/5 


BAYES NORM .91 


BAYES EP.1 .91 


BAYES E.25 . 9% 


BAYES TIDF .90 
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. 94 


- 94 


- 94 
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9% 


- 96 


~96 


96 


~95 
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96 
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EPSILON .1 


TABLE 7 


ERROR DISTRIBUTION 


572 573 574 5/75 


80 


.86 


86 


86 


~85 


-80 .78 


- Fim. 92 


908. 92 


790). 9X 


76 


- 9% 


- 94 


93 


oS 


EPSILON .25 


572 573 574 5/75 


66 


78 


- 78 


78 


77 


-61 .58 


-8¢ .87 


-8¢ .87 


.83 .86 


»82 .85 


WITH SCALES VARYING FROM .4 FOR #1 TO .8 FOR #5 


STUDENT T 1DF 


572 573 574 575 


“60 .57 


~78 .80 


-78 .80 


-81 .83 


-80 .82 


AVERAGE « 70 


WEIGH MEAN .91 


MEDIAN ~ 90 


VOTING 90 


BIWEIGHT 90 


MLE T 1DF .90 


MLE T 3DF .91 


Tr 


MEE GPS .1 .91 


MCE EPS.25 .91 


NOTE : 
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- 96 
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92 


a 
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1. PRIOR PROBABILITIES ARE ALL EQUAL 


2. FOR ALL ERROR, 


RECEIVER 2, 
RECEIVER 5 
3. THE SCALE OF THE SECOND TERM IN THE MIXED NORMAL WAS 10 TIMES 
THE FIRST SCALE 


4. ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
5. ALL ENTRIES WERE ROUNDED 


-6 FRO RECEIVER 3; 


41 


66 


66 


66 


- 66 


a7 


.69 


70 


~70 


OFF TO 2 DIGITS 


-61 .58 
78.79 
78 .78 
-77 .80 
81 - 84 
-81 .83 
-83 .85 


-83 .85 


SCALE PARAMETER WAS .¢% FOR RECEIVER 1, 
-7 FOR RECEIVER 4G, 


- 56 64 .61 
. 88 -71 .76 
89 71 .76 
87 -72 .78 
87 72 .78 
Be 60n.. 15 
53 63 .60 
5G. 63 .61 
81 -63 .76 
80 65 46 
81 -63 .73 
-86 67 .77 
~85 66 .77 
-88 -65 .76 
88 -65 .76 

-5 FOR 

AND .8 FOR 


-59 .56 


mo 419 


a5 .1& 


aio ./% 


ano .S2 


-78 .82 


/m 79 


-77 .80 


EE 


TABLE 8 
Pi TRUE SCALE GOES FROM .% TO .8, ASSUMED SCALE IS AVR 
ERROR DISTRIBUTION 
NORMAL EPSILON .1 EPSILON .25 STUDENT T 1DF 


ALGORITHM 5/2 5/73 5/74 5/5 572 573 5/74 5/5 5/72 573 574 5/5 5/72 5/73 57% 5/5 
BAYES NORM .91 .93 .94 .9% 79 .78 .75 .7% 65). Olas? 55 .63 .60 .58 .55 
BAYES EP.1 .90 .93 .94 .9% 68S 669 6 9hanaec 77 [82.65.07 70 .75 .7/7 fe 
BAYES —E.25 .90 .93 .9% .9%% ob -S89"-91 .92 77 Ses 185 1.87 70 .75 . 760 
BAYES T1DF .89 .92 .93 .9%% -85 .89 .90 .91 77 .82 .85 .87 .71 .78 .80 .82 
BAYES T3DF .90 .93 .9% .95 205 aoe. 9L 292 -77 .81 .84 .86 70 .77 .79 wee 


BAYES T10D .91 .93 .9% .95 .83 .88 .89 .90 ./2 .78 .80 .82 -68 .7% .76 <7 


AVE RAGE 91 .9353 .9%9 .94 79 .78 .76 «79 -66 .61 .57 .5% -63 .60 .58 s58 
WEIGH MEAN .91 .93 .9% .9% -79 .78 .76 .7% -66 .61 .57 .5% .63 .60 .58 J55 
MEDIAN 91 .96. 295 392 .79 .86 .88 .88 .66 .78 .79 .82 ~63 .75 .75 eee 
VOTING a9L «.9JOeeeo8. 9c o79° .6Gr 700" 707 66 .78 .78 .81 .63 .75 .25 “ee 
BIWEIGHT 91 .90 .95_.95 .79 .85 .88 .89 -66 .76 .79 F861 63 .735 .7/5 wae 
MEE@-USLDF .91>.92mr9S 2 oo -79 .<88).99 . 9a .66 .81 .82 .86 69 .77 af/9unee 
MEE T SDF .91 5955.94 5..9% 9 .66 7.90 <9. .66 .80 .81 .85 69 .76 .77 .80 
MLE.EPS .1 291 .95 274.7% «79 2.89 7072 .66 .81 .83 .86 .63 .75 .7655re 
MEE EPS-25 2.91 .92 .99 .% <M? .G9 .90-~92 66 162 388 .87 -63 .75 .7@ we 
NOTE: 
Ll. PRIOR PROBABILITIES ARE ALL EQUAL 
2. FOR ALL ERROR, SCALE PARAMETER WAS .% FOR RECEIVER 1, .5 FOR 
RECEIVER 2, .6 FRO RECEIVER 3, .7 FOR RECEIVER 4, AND .8 FOR 
REGEIVER 5 
3. THE SCALE OF THE SECOND TERM IN THE MIXED NORMAL WAS 10 TIMES 
THE FIRST SCALE 
G@. THE SCALE USED IN THE ALGORITHMS WAS THE AVERAGE OF THE SCALES 
FOR RECEIVERS IN USE 


5. ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
6. ALL ENTRIES WERE ROUNDED OFF TO 2 DIGITS 
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TABLE 9 
Pa PRIORS .15,.35.1».3,.15 RADARS ASSUMED .2 SCALE .5 


ERROR DISTRIBUTION 

NORMAL EPSILON .1 EPSLLON «co STUDENT T 1DF 
ALGORITHM 5/72 5/73 5/74 5/5 572 573 574 5/75 5/2 573 574 575 572 573 574 5/75 
BAYES NORM .69 .84 .92 .95 58 .67 .68 .69 45 .46 .45 .43 -39 .44 .46 .%6 
BAYES EP.1 .67 282 .91 .95 .60 .76 .85 .91 .50 .65 .76 .82 42 .55 .6¢4 .71 
BAYES —.25 .65 .80 .89 .94¢ soe) -/%e.G5 .90 -48 .63 .74 .81 -40 .53 .62 .70 
Beres FIDF .55 .75 .85 .91 -48 .69 .78 .86 39 .56 .68 .76 sae DUM? 469 
BAYES T3DF .60 .79 .88 .93 52 .72 .81 .88 41 .58 .70 .77 3% .51 .61 .71 


BAYES T10D .66 .82 .90 .95 -56 .73 .81 .88 -43 .57 .68 .75 oO) abe sabe 470 


AVERAGE 84 .92 .96 .98 Timo tr. 32 54 .51 .49 .47 -50 .50 .50 .50 
WEIGH MEAN .84¢ .92 .96 .98 lee Sel Salere 54 .51 .49 .47 -50 .50 .50 .50 
MEDIAN .8¢ .87 .93 .94¢ -71 .82 .88 .90 54 .72 .77 .82 -50 .69 .72 .79 
VOTING .8¢ .87 .89 .9G 71 .82 .85 .89 54 .72 .76 .80 -50 .69 .70 77 
BIWEIGHT .84 .87 .93 .95 7L wae .8e .92 -5¢ .70 .77 .8l1 50 .65 .71 .77 
MLE T 1DF .84 .89 .94 .96 71 .85 .90 .93 55 .76 .82 .87 5a .70 .768.82Z 
MLE T 3DF .84 .91 .95 .97 .72 .86 .90 .94 56 .75 .80 .86 “52 1.7/0 ./5'.81 
MEE EPS .1 .84 .91 .95 .97 «71 .87 .92 .95 .5¢ .77 .82 .89 .50 .68 .73 .79 
MLE EPS.25 .8¢4¢ .91 .95 .97 71 .87 .91 .95 54 .78 .83 .89 50 .68 .74 .80 
NOTE: 
le. RADAR 3 WAS ALWAYS SENDING 
Le THE ASSUMED PRIOR PROBABILITIES, BIRD 1 WAS .15, BIRD 2 WAS .3 
BIRD 3 WAS .1, BIRD @ WAS .3, BIRD 5 WAS .15 
ec FOR ALL ERROR, SCALE PARAMETER WAS .5 
G. THE SCALE OF THE SECOND TERM IN THE MIXED NORMAL WAS 10 TIMES 
THE FIRST SCALE 


5. ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
6. ALL ENTRIES WERE ROUNDED OFF TO 2 DIGITS 
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soo. dye Loe lo,.S ASSUMEO 


.2) SEATE.5 


ERROR OISTRIBUTION 
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P., PRIORS 
NORMAL EPSILON .1 
ALGORITHM 5/2 5/3 5/4 5/5 5/2 5/3 5/4 5/5 
BAYES NORM .80 .90 .94¢ .97. .68 .70 .72 
BAYES EP.1 .80 .89 .94¢ .96 .73 .83 .90 
BAYES E.25 .79 .88 .93 .96 .71 .82 .90 
BAYES T1OF .73 .86 .91 .95 .65 .80 .87 
BAYES T30F .77 .88 .93 .96 .68 .81 .89 
BAYES T100 .80 .89 .94 .96 .69 .80 .88 
AVERAGE 85 .92 .95 .97 .72 .72 .73 
WEIGH MEAN .85 .92 .95 .97  .72 .72 .73 
MEDIAN 85 .87 .93 .94¢ .72 .81 .88 
VOTING .85 .87 .88 .93 .72 .81 .85 
BINEIGHT .85 .87 .93 .95 .72 .81 .89 
MLE T 10F .85 .90 .94 .96 .72 .84 .90 
MLE T 3DF .85 .91 .95 .97  .72 .85 .91 
MLE EPS .1 .85 .91 .95 .97 .72 .86 .92 
MLE EPS.25 .85 .91 .95 .97 .72 .86 .91 


NOTE : 


nM 


RAOAR 3 WAS ALWAYS SENDING 


THE ASSUMED PRIOR PROBABILITIES, 


BIRD 3 WAS .1> 
FOR ALL ERROR, 


THE FIRST SCALE 
ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
ALL ENTRIES WERE ROUNOEO OFF TO 2 OIGITS 


BIRO @ WAS .15, 


71 750 
» 94 63 
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ae a5 
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92 ~55 
72 - 54 
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90 . 54 
89 - 54 
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SCALE PARAMETER WAS .5 
THE SCALE OF THE SECONO TERM IN THE MIXEO NORMAL WAS 10 TIMES 
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BIRD 1 WAS .3; 
BIRO 5 WAS .3 
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TABLE 11 


Pu PRIORS .05,.1,.15».3,.4¢ ASSUMED .2 SCALE .5 


ERROR DISTRIBUTION 

NORMAL EPSILON .1 EPSILON .25 STUDENT T 1DF 
ALGORITHM 5/72 5/73 5/4 5/5 B7f2 5435 57% 525 572° 573 5/¢ 545 572-575 5/74 5/5 
BAYES NORM .81 .91 .94 .97 70 “a5 33D al 52 .50 .48 .47 -48 .49 .49 .50 
BAYES EP.1 .81 .90 .94 .97 -/6 @B6 .91 94 65 .77 .84 .89 7560.67.75 .78 
BAYES E.25 .80 .89 .94 .96 -/5 .86 .90 .94¢ 65 “7/ «85 .89 56 .67 .73 .78 
BAYES THDF ./76 36/7 .92 .95 SicusS>S BBO 292 -63 .7¢4 .81 .88 .55 .68 .74 .80 
BAYES T3DF 179 .89 w9S ait ofS Gao 290 295 63 .75 .81 .87 56 .68 .74 .80 


BAYES T10D .81 .90 .9¢ .97 15 es 389 392 59 .71 .77 .84 29S 2.605 ./2 .17 


AVERAGE 84 .92 .95 .97 ~71 .7% «(67% «772 -5¢ .51 .49 .48 -50 .50 .49 .50 
WEIGH MEAN .8¢ .92 .95 .97 71 «7%: «.7G@ «272 -5G .51 .49 .498 -50 .50 .49 .50 
MEDIAN -84¢ .86 .93 .94¢ -/1 #61 286 290 -5¢ .73 .76 .84 -50 .69 .73 .79 
VOTING .84 .86 .88 .94¢ -/L «81 «85 .89 wO 275 «£5 .81 50 .69 .71 .76 
BIWEIGHT -8G¢ .86 .93 .95 77L Bl FSs8 292 7546 .71 7 82 -50 .66 .72 .77 
MLE T 1DF .8¢4¢ .89 .9¢ .96 -re .65 .90 .98 -54 .76 .81 .88 Sliver. 76 .82 
MLE T 3DF .84 .91 .95 .97 7/2 .85 .90 .94 55 .76 .80 .87 51.76 .75 .81 
MLE EPS .1 .8¢4 .91 .95 .97 -/1l .88 .91 .95 54 .78 .82 .90 -50 .68 .73 .79 
MLE EPS.25 .8¢ .91 .95 .97 ~-7lL #87 291 205 -5¢ .78 .83 .90 -50 .68 .73 .79 
NOTE: 
1. RADAR 3 WAS ALWAYS SENDING 
Za THE ASSUMED PRIOR PROBABILITIES, BIRD 1 WAS .O05, BIRD 2 WAS .1 
BIRD 3 WAS .15, BIRD @ WAS .3, BIRD 5 WAS .G¢ 
S FOR ALL ERROR, SCALE PARAMETER WAS .5 
G. THE SCALE OF THE SECOND TERM IN THE MIXED NORMAL WAS 10 TIMES 
THE FIRST SCALE 


5. ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
6. ALL ENTRIES WERE ROUNDED OFF TO 2 DIGITS 
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TABLE 12 


P, FOR AN END RADAR SIGNALING SCALE .5 
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NORMAL 


ALGORITHM 5/72 5/73 5/74 5/5 


ERROR DISTRIBUTION 
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-25 STUDENT JT 1DF 
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BAYES T3DF .92 .96 .97 .98 
BAYES T1DD .92 .96 .98 .99 
AVERAGE .92 .96 .98 .99 
WEIGH MEAN .92 .96 .98 .99 
MEDIAN .92 .93 .97 .97 
VOTING .92 .93 .97 .97 
BIWEIGHT .92 .94 .96 .97 
MLE T 1DF .92 .94 .97 .98 
MLE T 3DF .92 .95 .97 .98 
MLE EPS .1 .92 .96 .98 .99. 
MLE EPS.25 .92 .95 .97 .98 
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5. aie 
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RADAR 1 WAS ALWAYS SENDING 
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75 -75 .84 .87 .9D 


ASSUMED PRIOR PROBABILITIES WERE ALL EQUAL AT .2 
ALL ERROR, SCALE PARAMETER WAS .5 
SCALE OF THE SECOND TERM IN THE MIXED NORMAL WAS 1D TIMES 


FIRST SCALE 


ENTRIES ARE THE AVERAGE OF 10 RUNS OF 10D0 TRAILS EACH 


ENTRIES WERE 


ROUNDED OFF TO 2 DIGITS 
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TABLE 13 


Pa FOR AN INTERIOR POINT SCALE .5 ASSUMED PRIORS .2 


ERROR DISTRIBUTION 
NORMAL EPSILON .1 EPSILON .25 STUDENT T 1DF 
ALGORITHM 5/2 5/3 5/4 5/5 5/2 5/3 5/4 5/5 5/2 5/3 5/4 5/5 5/2 5/3 5/4 5/5 
BAYES NORM .84 .91 .96 .97. .71 .72 .72 .72 .54¢ .51 .48 .48 .49 .50 .50 .49 
BRvESuEne) .Squ,9195u897 © .78 .87°.92995 .70 .79 .85 .90 .61 .69 .75 .79 


BaAves —E.25 .83 2908.95 .97 wre .8/ .92 294 70 .79 .86 .90 -61 .70 .75 .80 


BAYES TIDF .82 .89 .99 .96 7/7 .85 .91 .99 -69 .78 .8¢ .89 "62 2.7L .70o .82 
BAYES T3DF .8¢4¢ .90 .95 .97 -/7 .86 .91 .99 -68 .77 .8¢4 .89 CULE T? .O2 
BAYES T10D .8¢4¢ .91 .95 .98 [75 .8635.609 2.935 ,Osee12 «795265 56 .68 .73 .79 
AVERAGE 89 .91 .96 .97 ml «iewere w/e -5¢ .51 .99 .98 -4¢9 .50 .50 .50 
WEIGH MEAN .8¢ .91 .96 .97 Sey eee? Come) & > -sutae -54¢ .51 .49 .48 -49 .50 .50 .50 
MEDIAN -89 .86 .93 .99 s7L sSh .88 .90 5¢ .73 .77 .83 -49 .69 .72 .79 
VOTING -8¢ .86 .89 .94¢ 7/1 .81 .85 .89 -5¢ .73 .76 .8l1 49 .69 .71 .77 
BIWEIGHT 89 .87 .93 .95 7/1 .81 .88 .91 -5¢ .70 .77 .82 -49 .65 .72 .76 
MLE T 1DF .84¢ .89 .99 .96 -71 .84 .90 .93 Beem sO ol .Or leer? O 1.76 .81 
MLE T 3DF .8¢ .90 .95 .97 -/l .85 .90 .93 Son .75 .GUn.or -51 .70 .74 .8l1 
MLE EPS .1 .84 .91 .95 .97 za -86 .92 .95 59 .7/ .82 «89 -949 .68 .73 .78 
MLE EPS.25 .84 .90 .95 .97 -/l .86 .91 .95 5% .77 .82 .89 -49 .69 .73 .79 
NOTE: 

1’ RADAR 3 WAS ALWAYS SENDING 

2. THE ASSUMED PRIOR PROBABILITIES WERE ALL EQUAL AT .2 

Se FOR ALL ERROR, SCALE PARAMETER WAS .5 

G. THE SCALE OF THE SECOND TERM IN THE MIXED NORMAL WAS 10 TIMES 

THE FIRST SCALE 
5. ALL ENTRIES ARE THE AVERAGE OF 10 RUNS OF 1000 TRAILS EACH 
6. ALL ENTRIES WERE ROUNDED OFF TO 2 DIGITS 
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2. MULTIPLE PARAMETER CASE 

The following tables are the results for the two parameter simulations. Each 
table represents a particular case of interest. The entries are the percentage of correct 
identifications for 2000 trials. 

The tables contain several different conditions each. Across the top of each table 
we see that the first subdivision is rho. This represents the value of rho used to 
compute the error or “true” rho. Rho ranged from -.5 to .5. Directly below the rho(p) 
are the number of radars/ number of observers used in to both generate and analyze 
the data. The left hand column is again broken down into two parts. The large 
subdivision is the value of epsilon used to produce the errors.’ This ranged from 0 (the 
bivariate normal) to .25. The next subdivision is the particular algorithm used to 
identify the radar and it’s associated epsilon or degrees-of-freedom component. 
| In order to understand the abilities of the particular algorithms we examined a 
large number of situations. We looked at two patterns, box and diagonal, under 
various condition of changing true and assumed( what was used by the algorithms) rho 
and sigmas. 

The tables below contain the results for the probability of hit. They are divided 
into two large groups, the first is for the diagonal pattern and the second the box 
pattern. Each pattern based group 1s further divided by the scale used to generate the 
error, ‘true scale”, The combinations of scale parameters for the x and y components 
were (.5,.5), (.3,.7), and (.9,.9}). With in each set of true scales there are difleremmsemies 
used for the algorithms, “assumed scale”. The true scale set of (.3,.7) were (.5,.5) and 


(.3,.7). Finally, all the combinations were evaluated for rho of -.5, 0 and .5. 
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TABLE 14 
Pi, FOR A SCALE 0.9, ASSUMED RHO = 0.5, LINE PATTERN 
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TABLE 15 
Pi FOR A SCALE 0.5, ASSUMED RHO 0.0, LINE PATTERN 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 16 
Pi, FOR A SCALE 0.5, ASSUMED RHO OF -.5, LINE PATTERN 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 1/7 
Py SCALE .5, ASSUMED RHO OF .5, LINE PATTERN, RADAR 3 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 18 
Py SCALE 0.5, ASSUMED RHO OF 0, LINE PATTERN, RADAR 3 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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Pi, TRUE AND ASSUMED SCALES OF .3 AND .7, RHO .5, LINE 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 


RHO = -.5 


RHO = 0 


=o 


RHO = 


EPSILON=0 
A 
M 


i) 
Ne leat tenet 
i) 


Seta 
in) 

COD DA 
SAA AAADO 
i> e e e e 
QyDI~ 0000000000 
LAAAHAAAHGA 
if) 


ee @ @©@ j@  @ @ # 


ine 

ener ete te tee 
il> e 

+O? 

LDA ase 
i> ° 

cn 0 
eee 
it> > 

Ay 000001000000 
SHAHNAAAAD 
iw) 


LN 
SANs 
im: 

ph 
SNe tet 
im: 

ot 
SOA 
im: 
QyTI1ACOCIDA0O 
SAAAAAAA 
iw> 


e¢« oe# e«® @e®  @®  @e  @ 


a a — 


EPSTLON=. 1 


in Ost 
SDM eaedeset 
it> e e 


ph st DO 
SS DOMAAD Ae 
it> © e e 


ODES DACOAAD 
SDONAAADA 


qh mcolrnwor~ 
S_WUONHAAAAA 
iw) 


pn COst 
— OrcOe sete tert 
i> e oe 


fFOLNAAAVHA 
ST NOAAAAA 


on Nst a000000 
VC DAOAAAAA 


QySt Lon ston 
SOON AAAA 
ir) 


LM LOvo 

SD MA rarstedtes 
i> s e 

pti NA 

NS DOAHN Arar 
iu) s & s © 


HAWWAACOAD 
SCOCID MAA 


QITMWOUWLNWW 
LUO CNAAHA 


>. ee: @®  @ oe #@  @ 


a 


=.25 


EPSILON 


LQ OOM MOADA 
SS DONAAAGA 


<-0)000000\9 0000 
SDWDOANAAAA 
Lm 


MOALALNMLOLN 
COMMMMATAGY 


QIN AO etes 
SIS HAMDCOMA 
ir) 


LATION MNOOMA 
SS DODAAAA 
if) 


<p} LOP~CO0OM 0000 
LS DWAAAHAA 


LA COWLNMWiLO 
SCOONDMAMMO? 


QMAtOONr4 
LOM AACODA 
Fe) 


LQNTHIOANANAD 
SDI AAAHAA 
if) 


>. @ @  @  8@  @ @ 


jpuONWNAUOAA 
ST DOAAAAA 


MIAWWst Www 
SLOW AAAAA 
in) 


qy = AHOOne 
LWMAACOAA 


=< —F_F—_UL 


CAUCHY 


LA NOOCO0Om mm 
SC HDWOAMNAAHA 


pO WWW 
2 COM AAAAA 
un 


MNOSFMNMA 
COM DAAAMA 


QyDRLOLN sHLOLN 
~~ Wr.00.000000 00 
Ln 


NONnrnmwonrm 
SDI AAMAAG 
in) 


St NLL st tt 
LOM AAAAGA 


OAINNNOOO 
COM AAMAAG 
uw 


QOOttnus¢m 
~SWor.00 00000000 


LOh= NESEY OO 
S_COPAAAHGA 
if) 


>. e@ @  @  @  @# ©@ 


pONMMAHt tH 
SCOP AHAAOIO? 
iw) 


MQOANNHOO 
SSIS AMAMAAG 


QyTtNWOWWMLMWY) 
Svwor.000c0000000 


WO => A.a0- 
=>—ZFFFE UU 


D 
N 
NS OF 200 TRIALS EACH 


6] 


TABLE 27 


Pi) TRUE AND ASSUMED SCALES OF .3 AND .7, RHO O,LINE 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 


RHO = -.5 


RHO = 0 


2 


RHO = 


=0 


EPSILON 


nS” 
SOA 


Lo ® 


Yat 
uM 

m0 
SOA 
Lo e 

ay 
SOA 
Lo e 


Ln 
SOA 
Lm e 
<0? 
SOA 
Lo ° 
on 
SOAs 
Lo ° 


QJLODAAAAO 
LAMDAAAAGD 


LOCO 
SOA 
Lo ° 
<p h™ 
SOA 
“tom e 


OMIAAHAADHO? 
NK AAAAHA DIO? 


EPSILON 


incor 
SN DNO A ases 
LO e ° 
hh 
NOAH 
Lm g e 


atT DAD 
SNOADHA 
Lm e e e e 


yh =QO000Om C000 
SCOCDAAADA 


LO h=sO 

SN HDOO MAA ete Ate 
Lm e e 

«pz LOO OY 

SN DOAANAA 
Lm e e e 


ONS ADAH 
SDCODAAAAO 


qs Orr vr r 
S_COCTDDA AAO 
un 


ineco 
DS DNO AAA 
Lo e e 


op FLOOD? 
SN HADOAAD AA 


cr DLO COCO 0000 
SE OCODAAAA 
LO 


QIN Wwlnww 
WS COCTDM ATA TAG 


=.2 


EPSILON 


LNTON Or 
SOS AANA 
L- e e * 


<p DUADIDIOND 
VON AAAACH 


ey) DYNICOCOM™ 0000 
SOM DAAAGA 
LO 


QyNWtMALOLO 
SRAM AAG 


LOANMAAAAA 
SOI DAAAOY 
Lm 


<p OODDMAO 
DT AORHDAAAG 


7 ONIPSLOLORS ES 
OMS AAAAON 
LO 


QOtNNON™M 
SSIS DAA AO 
Lo 


LN TINAAOAA 
SON AAAAGA 


<j M0000W0000 
SCOP HAADHDIO? 


AMINWOWAFWW 
NCOP AAAAO 


CAUCHY 


LN NWODAAHAA 
SOI DAAAG 


<p LOTVCOCOMSM 
TL COMAAAAG 


or) OP WWOLOLONWO 
COM DATA 
LO 


QQKNOANAOoO 
LAR AAAAG 


Ln PLM COCOM 
SDAIN DAAAAG 


<p LOLNWOWOWOLOLD 
SOM AAAAA 


MOUWTFAMMM 
SOM AAAGAGA 
Lo 


qQyCOWW CICOF~ COM 
SLOP 00.00.0000 CO 
uM 


LN MWOCMMmWOWWw 
SK CORAAAAG 


g$NOTIMMM| 
DCO DAAAG 
LO 


mhTAOONOCO 
SPY PF. 090000000 
uM 


VMVNOtmMaAa 
SLOP 0000000000 


=. <_t—F-F- LL 


NOTE: 


D 
N 
NS OF 200 TRIALS EACH 


62 


TABLE 28 
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ele 29 
Py TRUE SCALE .9, ASSUMED SCALE .5 , ASSUMED RHO 0.5 
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TABLE 30 
Py TRUE SCALE .9, ASSUMED SCALE .5, ASSUMED RHO 0.0 
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Py TRUE SCALE .9, ASSUMED SCALE .5, ASSUMED. RHO .5 
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TABLE 53 
Pi, TRUE AND ASSUMED SCALE .9, ASSUMED RHO 0.95 
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TABLE 33 
Py TRUE AND ASSUMED SCALE .9, ASSUMED RHO 0.0 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 34 
Pi TRUE AND ASSUMED SCALES .9, ASSUMED RHO .5 
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WABLE 35 
Pry SCALE .3, ASSUMED RHO 0.5, BOX PATTERN 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 36 
Pi, SCALE .5, ASSUMED RHO 0.0, BOX PATTERN 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 37 
Py SCALE .5, ASSUMED RHO -.5, BOX PATTERN 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 38 
Pi, SCALE .5, ASSUMED RHO 0.5, RADAR 3, BOX PATTERN 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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Pi SCALE .5, ASSUMED RHO 0.0, RADAR 3, BOX PATTERN 
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TABLE 40 
Pi SCALE .5, ASSUMED RHO .5, RADAR 3, BOX PATTERN 
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TABLE 41 
Pi SCALE .5, ASSUMED RHO 0.5, RADAR 5, BOX PATTERN 
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TABLE 42 
Pi, SCALE .5, ASSUMED RHO 0.0, RADAR 5, BOX PATTERN 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 43 
Pi SCALE .5, ASSUMED RHO .5, RADAR 5, BOX PATTERN 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 44 
Pi, TRUE SCALE .3 AND .7, ASSUMED .5, RHO .5, BOX 
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Pu TRUE SCALE .3 AND .7, ASSUMED.5, RHO .0), BOX 
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TABLE 46 
Py) TRUE SCALE .3 AND .7, ASSUMED .5, RHO ~-.5, BOX 
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TABLE 47 


Pi, TRUE AND ASSUMED SCALES .3 AND .7, RHO .5, BOX 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 48 


Pi, TRUE AND ASSUMED SCALES .3 AND .7, RHO .0, BOX 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 
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TABLE 49 


Pi, TRUE AND ASSUMED SCALES .3 AND .7, RHO 


=) OA 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 


RHO = -.5 


RHO = 0 


5 


RHO = 





=0 


EPSILON 





uO? 


SDAA 


uo e 
tO 
SOMA 
wo e 

oO 

SOS ane 
wo e 


@yst000000 000000 
S—AHNAHOAO 


uO? 
SOA Maat 
um ° 


pfAOAAAAOO 
TL DAAAGAGOG 


cn OMCOAAACD 
SAAAHNAAAOY 


QypTowWMwoiwuwnn 
NS AADAANIHAOIO 
um 


LQNCODDADAAANO 
SK DHAAAAAG 
Te) 


=: *§ *  @  @  @  @ 


<p VAIMADACOO0O 
NS AHDNIANODIOY 


cn LOR PP COP oO 
SAAD OIOIOIO?Y 
um 


Quarts stan 
SANMAAG OO 


=< t—F-F-LULU 


rc 
I 
—_ 
© 
=) 
I 
~Y 
CO 
Li) 


wor 

VS DOMAS4aeste4 
wom e e 

< OM 

SN DO MAMAS 
wm = e 

on TIF OONCODO) 
SHNOAAAAO? 
Lm 


Qype OWMWO 
SDOAHADAAOG 


LH ONDAAAMA 
SHDOAAAAG 
un 


qpiOCIOADAIMNdy 
SHDOAAAAO? 
i @) 


nw stern cors 
SHDODHAIAAOD? 


QysPLomma§iman 
WS COCOMGMAAO 
Te) 


LOh= MOON CO 
SDOHNAIAAA 
“Wm 


<p LOQIPSPSPSP WO 
SHAOIHAAAAO 


MVOstst tt 
CHNAOAAAAG 
uM 


QM s4sHOON 
~~ OCMNAAAICN 
“uM 


= @ @ #$j@  @  @ & 


=< -—-t—F-LUIL 


©) 


=.2 


EPSILON 


LON StF DDOMO 
SS DWDAAG OA 
un 


jpOMOMOWDNO 
S_WDWODHAAGIOY 


MOWWMMWW 
COWNAADG 


~_ LOM HCOCOIOY 


L727 COCOM.0000 
KN OODAAODOY 
ut) 


<¢ COLMPLOrMmrn 
SOW DAAGA OY 
ut) 


omhstston9 st 
S_COLWDMA AAO 


Qeor water 
SPW O00 00.00 
Lm 


LN COSI Ovrom 
SHO AADOVOY 
i @) 


—jf~OMNUOMNLO 
SOW DAAAOY 
i @) 


on FuoOoOo0nd0rnoy 
~_COW DIDICN COCO 
um 


QyComwostorst st 
SOW COCO COCO0O 
Lo  e@ ee e# oe @® @ 


= —tF-F-F-LUILY 


CAUCHY 


LO MUM WWWWw 
~COMAAAAO 
Te) 


pPONOMAMO 
S_COMAAADVOY 


QyCONILOLO tO +t 
~SLOr-000000.0000 
un) 


InNowunuststt 
CS NWAAAAG 


SOWDAAAG 
unm 


on SOO GOO 
~0010.0000000060 
uw) 


~_Sowaaoorm 
“um 


NOMTHTMMNM 
CHDUOAAAAG 
Ln 


plot AMTAOOO?) 
SS COWOGC00N00 
un) 


es us “eo Ss “@ "6" “4 


OQWINMWWOTFM OY 
~00.0.00.00.0000.00 
im 


QHOMDAAMmWrwo 
SOON NRE 


== t-t-F-Lutu 


NOT E: 


WZ =i Lu 
al) ee ae 
Sr FR F 
t—- e414 
ox eeer ze 
Or OO re 
EO CE 
au. «= J OE 
Lu <IM™ Lu 
™ ~ WY 
Lu Lum = 
i 22 ae 
=z OF bh 
<< © ee 
Lu a €©6C 
——_~ <i t 
<C . ©) 
MO oct + 
WIS Llu le 
=zt Me tL 
Lu Diu LL 
Com = €) 
li) ti<t Cc) 
Or —-Iom 
FLU <C<f 
ae Qa. 2 
Or NY © 
LI t- 1 
M<rl LUN FF 
JO eae < 
tr Ss 
Lut— i 
—IY) @@e oc 
<r Ow me 
Cre Ur oO . 
“ALL i Comm 
LL e 
LLU oe COLL I 
ae 
H--+-F—-<tu FF GS 
uj <C 
Lica ~~ Ere 
Ooazwo OF z=. 
LoL wo Mart 
lj <> &- lus 
IMAG Dak 
el ZN —J}-—t— 





84 


TABLE 50 
Pi TRUE SCALE .9, ASSUMED SCALE .5 , RHO 0.5, BOK 
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RHO = -.5 


RHO = 0 


<6 


RHO = 


=0 


EPSILON 


INOTAWWONeH 
~_0007000000.0000 


<p COCONIt 0000 
~ cococococomnm. 
ut) 


MOMOMAME) 
~~ 0Omm. comm 


QyT~manndtst 
SNR WOWrNWOW 
un) 


NONRAATA 
~00.G10000070000 


piOCMMcaand 
~_00.00 COCOCOr ™ 


MONCOMAAILN St 
~oomoocomm 
ir) 


>. @®° @  @  @ oe 


Quy taaMrwo 
SE ONE.WOW 


uM ._ e® se @® oe @® 


in OMOnTANAIOCOm 
SN COMDIG HON 0O 


jpfOTWOOWOst 
00.0100 09.5100 CO 
im 


MOrNtWOHO 
~P000000000000 


>. © © $@® ese @e@®  @ 


= FF} LUI) 


INWOtONTHO 
~_0cOMm.00.00000000 


pNHAONTOWM 
SCO 0000 COM 


MHODMN OOM 
Sh Orem 
um 


Qh e CWO HAOIAI 
N_WWOWWWDWWWO 


Lp TON BOM 
~_00\000 00010000 
ir) 


= @© @  @  @  @  ©@ 


~_COrP~00 CO. COr 


nHtatwouno 
CM Oonnnne 
Ln 


qQyCOP-CODMMmWN 
SWWOrEWW 


Int AWwnOWMSt 
~_00™.00.0000.00.00 
it) 


= e¢& @©&  @  @  @  @ 


gNOtonton 
~_CO00. 00000000 
if) 


MMOORTAOAIMWO 
Sr CO COCO 


= j® @  @  @  @  ©«@ 


aa 


=. 


EPSILON 


NW DHDAANE 
Streococomer. 
in 


> ses &©® &® @® ae 


pHIOTONHO 
SLO NE 


MONDO MWOwW 
COLNE OW 
in 


OO. Ooo Pies OW: 


QHOmMANst 00 
SNM WOWUOMLY) 


Ip OMaAMMcaw 
~™st cococomm 
it) 


> ¢ @ @® @e®  @ @ 


prinwododtas 
Stor re COPS 


MOFAOnMMOrWN 
SOLON WOW 
un 


QJNMMNWHAO 
CUM WOWWWW 


IQNHNwotwouno 
~_M™st+0000000000 
ifm) 


= e&¢& © @  @  @ 


pfOIMMNAIMO® 
SPS 00 00. 00. CO ™ 
“im 


oh MMF COW) 
SOW NSS 


== FFF —- Lat 


CAUCHY 


instr oonm tt 
Steerer 
it 


jpfAOHAIMCINOAO? 
SWOT PES ESO 


mMmorcrmiust 
TWMWWWWW 
in 


= j@® @  @  @  @ 


qyhNaooonrn 
LUN tWOWWNIN 


INNOr OOM 
So mermereserer 
“Lo 


so ce e®# @&@®  @  @  @ 


gph HAOICIAICQCO 
NWO TTT IS WOW 
it) 


MNDOOMOM 
OWE EWOWW 
un) 


qyhorOooor. 
SIN TWOWWMNLN 


NMANMNAO 
~™in 0000000000 


ph NE COW tH 
SOW ESSE 
“tm 


MOMMMNNOO 
SOU NAN 


= &¢§ @  @  @e@  @  @ 


LUO > = = 0.0. 
=< F—-F-F-U 


NOTE: 


o> Lu 

— <I as 

<LY -— 

CYL a 

ES en si 

cock 

or Ce 

Lisle Lu 

cr WY 

Li<t — 

os 

-— -— 

© al 

Ltsrw— Lry 

t—f-— — 

<C<t > = 

orc) i Co) 

Lijit Lit << 

pe Lil Li 

Li jm LLU 

CS © WY 

—_ Rite 2a) 

Oul <— 

Ro © 

MH Cra 

- <LO - 

LLC yD et 

OZ ke 2) 

De<t = @- 
ON 

LISLIUSOCLLsLey 

IS 4eroc LL 

It <o © 

O~us Oo 

MW HfMNMO WY 

O< rly Nom ZZ 

Lts3k- —JLL — 

ae cr 

H-C3OLLIC IFN) 

Wed = <t 6} 

Um us re 


>< on. <L 

Oo & 
LILO LULL COS><Lu 
ge oe ae a a 
rere oot 


a Nat ww 


SS 


TABLE 51 
Py PRUE SCALE .o) ASSUMED SCALE. 3, hn@eO: 0 EUR 


CORRELATION COEFFICIENT USED TO GENERATE THE ERROR 


RHO = ~-.5 


RHO = 0 


= .5 


RHO 


=0 


EPSTLON 


LN WMAaaiocom 
S_COMICOMGICOCO 


j+~WDOWWOWS 
COM CONCOM0000 
im) 


e¢- @e@ @# @  *  @® @ 


7 AMWMOO 
~~ COCOCOCOCOM 


ay hYotdte AO 
SPH PY Preps pepe 
Lm 


INOMMAOWMWM 
~_COMOCNGIM 000 


qpfOoAMWOdre 
~~ COCO COCO CO COCO 
Lo 


e¢-06UM8RmUr OU SS 


MOS DHAMCW! 
Sr cor cocorm.P~ 
o> 


Qyh Yat wonm 
SIS PSPS SISO 


INOIMMNOAQICW 
~_COMCOMM COCO 
ut) 


j~DWOOrrASFMY 
~CO.TD1ICON COCO COCO 
LM 


= e¢ @  @  @  @&@  ¢ 


MAWAIMNOODN 
~_Po0 co. 0co.co0com 
“Lo 


QyAMMWcdac 
Se 
im 


=z —F—- FLU La 


QOS NOMWKN 
~_COM 0000010000 


geNOtRmMe 
~ cOm0000000000 
Lo 


MONHAONCW 
~rmrcocomem 
o> 


QyP =P ONMOO) 
SWUWRNE ERO 


NtaAwOMoOMs 
~_0OM00000700 00 


S comcoco000000 
“o> 


mtawoMoomnst 
SSIS Roomre 
im 


SWOOPS PSPS 00 
i> 


IN StrWAMING? 
~_COM0000000000 


S cor~0000000000 
if> 


mM AWOO-wMN 
SIM Pc000mrn. 


QyHMARANS HO 
SWORE 


= 2 


EPSTLON 


SPN C000 CO.00.00 
tf) ee e«# @  @  @  @ @ 
ptt 
SLOP. 
Lm 


OAMFATPTMNOMAI 
SWOPE 


qyOwrvworooun st 
SLNINWWWWW 
1S @) 


NWI NSNNO 
~_Pt+0000000000 
o> e — *° oe 
gi toortmwn 
~~ wor. co00r 


MOMNUTNAO 
SOLON AIRES 
Ce) 


QyMLOLOP~ COL) 
WN MMWWWWOwo 


LNNCOM SLOT 
~~ t coCoCOCOM 


Swe COOOrr 
im> 


oir NFWOUWTM 
~LOLO RAN 


CAUCHY 


LN SATAN 
~_ Im oocoorrm 


tft? e¢- ee: @  @  @  @  « 
qGOMONMuMuNt 
SOLON 
i> 


MOM ANAICIDIO 
~_LOLOEPY PO’ 


QGQNnnrdsdtoc 
UNNWWWWW 
“o> 


NNTIONnOMMN 
—™Lormroomrm 


gh Oss 
SOLON 


mONTAOOOWD 
~_LOLOESP PF Ow 
LO 


QyQOONINIRDNO 
SNMWWNWUMNW 
if) 


LNMASHONA 
~in.coco00m me 


gl-COWOUlOMOMEY) 
~_ WOW ESP 
1S @) 


mNFAHAOONW 
SLO LOPES PS WO 


qygTttt3AIAOO 
WN MMWW Www 


NOWE: 


ea Li 
=< ele, 
IY -— 
crcl Pal 
)=I- -— 
Oc 
oc Ee 
LrjLiL Lt 
ox Y 
Lu<t > 
=~ 
— 
© — 
Lijn— Lu 
-—-— — 
<<t © alg 
ocd) 4 ©) 
Lrj~ Lie <U 
wee Li Liu 
Lrj-~ Lt 
Cor a) WY 
al <)> ah 
OvrI mM < 
tC tH 
— Zaz co 
OO MZzO a 
LrjlzJ <Cr4 
NTO kt 7) 
— rz ) 
<[ AJ 
LrILi Lt 
is Luo? Li 
I<t .rcm cS 
Os to 
Mm lI cw WY 
Mcaciu = 
lk—<l LOO] 
ae <C Cc 
H—COLLUCOFE—WY 
Wen xT CC CE 
Uuiwm<at US 


a ee ealnd 


a4 NM ww 


86 


TABLE 52 
Pi TRUE SCALE .9, ASSUMED SCALE .5, RHO .5, BOX 
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Py, TRUE AND ASSUMED SCALE .9, ASSUMED RHO 0.5, BOX 
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TABLE 54 
Pi, TRUE AND ASSUMED SCALE .9, ASSUMED RHO 0.0, BOX 
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Se 55 
PH. TRUE AND ASSUMED SCALES .9, ASSUMED RHO .5, BOX 
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APPENDIX B 
SIMULATION CODE FOR SINGLE PARAMETER CASE 
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4 Met AG ent (15), ITRUE) 


DRAWS THE RANDOM NUMBERS AND COMPUTES THE 
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1,NRUN 
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Go 


250) ,REPMN(4,4,15), EPS(10) 
.THE FOLLOWING DO LOOP DOES THE RUNS 


STON DSEED 


DO 250 NWATCHs 


CALL DRAW(DSEED,TEPS) 
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.INTIATIALIZE COUNTERS FOR THE NUMBER OF CORRECT ANSWERS 
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BO J00R TREES 


10 
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CONTINUE 


CONTINUE 
.MEAN NUMBER OF CORRECT CHOICES FOR 1000 RUNS BY METHOD A 


WRITE(6, 300 
300 EORMAT CG ates 


301 EORMAT eI BAYES 


iS 
.NOW WE COMPUTE THE MEAN FOR EACH METHOD 


OUTPUT FOLLOWS 


ou 
100 CONTINUE 


250 CONTINUE 
350 CONTINUE 
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(M 


( 
WRITE(6,320) ((REPMN(M,N,14),N 


307 F 


,4( 4083. 2 UD) 


320 FORMAT (' WEIGH MEAN 


age" 


WRITE(6 
FORMAT -( 
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32] 
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WRITE 
310 FORMAT 
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* SUBROUTINE INPUT 


KKKKKKK KKK KKK 


Se CONTAINS ALL THE COMMANDS FOR INPUTING THE NECESSARY 


by2( 10) 


4c 


. ENTER THE NUMBER OF BIRDS 


NBIRDS=5 


Seienestme VALUES FOR EPSILON 


EPS(1)=0 
EPS(2)=. 1 
Emit )=. 25 
EPS(4)=4 


. INPUT THE MEANS 


TMEAN(1)= 1 
TMEAN( 23= 2 
TMEAN( 3)= 3 
TMEAN( 4 )= 4 
TMEAN(5 = 5 

INPUT THE A PRIOR PROBABILITIES FOR EACH BIRD 
TPRIOR(1)=. 2 
TPRIOR( 2)=. 2 
TPRIOR( 3)=. 2 
TPRIOR( 4)=. 2 
TPRIOR(5)=. 2 
er INPUT THE PARAMETERS FOR THE SELECTED ERROR DISTRIBUTION 
TSTDV1(1)=. 4 
TSTDV2(1)=" 4 
TOF(1)=1 
TSTDVA(2 =.5 
TSTDV2(2)= 5 
TDF(2)=1 
TSTDVA¢ 3 =. § 
TSTDV2(3)= 6 
TOF(3)=1 
TSIDVA(4)=. 7 
TSTOV2(4)= 7 
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.. ANALYSIS PARAMETERS 
ENTER THE NUMBER OF BIRDS BELIEVED TO BE IN AREA 
IBIRD=5 


ENTER THE ASSUMED MEANS 
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X(I) 
CONTINUE 


80 


ELSE 
DO 100 I = 1,NWATCH 
COMPUTES ERROR BASED ON EPSILON CONTAMINATION OR NORMAL 


RAN2 = GGUBFS(DSEED) 
CALL GGNML ra ED, 1, RAN3) 
IF ( RAN2 .GT. TEPS) THEN 
cpelth = RAN3*TSTDV1(I) 
DELTA = RAN3*TSTDV2(T) 
END IF 


X(I)=ZMEAN + DELTA 
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(DI+1)*LOG(1+(((X(1)-AMEAN(J) )/ASTDV1(1))**2)/DI) 
QR} Y4(3) 


ROBA 


NWATCH 
ITRUE) CNTA1= CNTA1+ 1 


1, IBIRD 
_ 


Q 


SUM=SUM 


CONTINUE 
XA . EQ. 


0 

0 

J 
THIS IF STATEMENT PREVENTs TO GREAT A ARGUMENT FOR EXP 


ao 

TA 

100 

SUM 

DO 50 [= 1 


SUBROUTINE TBAYES(DI,CNTAL) 
M 
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REAL UA(10) , PROBA 


SU 
TE 
DO 
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100 CONTINUE 


. 5-MULT DEE 


PROBABILITY 
OUT Tit 


STERIOR PROBABILITIES 
5 WRITE Pp 
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SUM2 = SUM2 + (((XC(I)-AMEAN(J))/ASTDVICI))**2) 
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CONTINUE 


eR EO ee AO EL 


DO 50 I = 1,NWATCH 


SUBROUTINE NORMBY(COUNTC) 
SUMZ = 0 


DEFINITIONS OF LOCAL VARIABLES - 
REAL UB(10),PROBB 


C 
e 


COOODOO00O 


50 


100 


1, IBIRD 


Doe ZU rd 


CHECK FOR NUMBER THAT WOULD BE TQ LARGE. FOR EXP 
TESTB = 0 


100 CONTINUE 


C= 


END IF 
200 CONTINUE 


TF_( MAXB .EQ. ITRUE) COUNTC = COUNTC + 1 
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443) )-¢a/(astovet1y)**2) ) 
2(1))+(. S*(DIFF/ASTOVZ( 1) ee 


.EQ. 0 ) WHEN 
100000 
COUNT+1 
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) THEN 


2 


LOR 
U(d 
PU(J)=LOG(APRIOR(J))+ PU(J) 


GT. RMAX 
PU(J) 
ITRUE) COUNT 


IT = 1,NWATCH 
BuSE 


MAX=J 
END IF 


50 CONTINUE 


= 


DIFF=ABS(X(1)-AMEAN(J) ) 


1, 1BIRD 

3 

ia ( ae 

END. 1G 
Pa 


PU(J 
DO? 
MAX . EQ. 
U 


TE 
REY RN . 
END 


SUBROUTINE EPSBAY( COUNT , EPS) 
CONTINUE 


REAL PU(20) 
RMAX=-10000 


pO 30" 
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* SUBROUTINE MEAN * 
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DEFINITIONS OF LOCAL VARIABLES 
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THE OBSERVATIONS 
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COMPAR(AMEAN, IBIRD, ITRUE,AMEDIN , COUNTE ) 
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RMEDIN(X,NWATCH) 
pCeereomOeeeemm: MEDIAN WITH IN .5 OF TRUE VALUE IF IT IS ADD 1 


= 
== 
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SUBROUTINE MEAN(COUNTOD,COUNTE,COUNTF) 
URN 


COUNTE 


COMPUTE THE AVERAGE AND THE WEIGHTED MEAN OF THE OBSERVATIONS 
AMEDIN 


SU) CONIINUE 
. COMPUTE MEDIAN 


C>. 
Ca: 
oe 


TON NO X(I) .GE. UCN)-.5) THR 


COU 


N) .£Q. MAX1) THEN 


) 


THE MAXIMUM NUMBER OF HITS, IF 2 OR MORE 
AX1) THEN 


bolnT@dy#i 


"( 


2 
LE. 
GT 
OUNT(N) 


i 


X,U,NWATCH, IBIRD,COUNTG, ITRUE) 
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IBIRD-1 


I 


? 


is 


( 1 eels 
cpuNTCY) 
in 


1,NWATCH 
C 
G 
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N=2 
PF 


COUNT(1) 


IP, {l 

] 

=0 

END IF 
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b ITRUE, AMED ,COUNTG) 


ITRUE ) COUNTG 
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Ease 


COUNTG+] 
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IF ( MAX . EQ. 


END Tr 
URN 
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~ SUBROUTINE BIWGHT 
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THIS SUBROUTINE IDENTIFYIES THE BIRD BY THE METHOD_OF BIWEIGHTS 
PESCRIEED IN BERERENCE ©. WRAL 1S GENERATED THAT IS 
THEN COMPARED TO THE GIVEN MEANS OF THE BIRDS 


: DEFINITIONS OF LOCAL VARIABLES 

C ALSTOP. - COUNTS ITERATIONS AND STOPS IF EXCEEDS 1000 

‘ DIFFCI) - RETAINS DIFFERENCES BETWEEN OBSERVATIONS AND UHAT 
C SUM - USED TO_SUM NUMERATOR 

C SUM2 - SUMS THE DENOMINATOR 

C UHAT - THE CURRENT ESTIMATE OF THE TRUE MEAN 

C Ul - RETAINS MEDIAN 

C U2 - THE LATEST ESTIMATE OF THE TRUE MEAN 

C W “THe WEMGHT 


SUBROUTINE BIWGHT(BIWAIT, IFLAG) 
COMMON AMEAN( 10) APRIOR( 10), ASTOVIC 
COMMON IBIRD, ITRUE,NBIRDS NWA 
COMMON TDF(10),TMEAN(10),TPRIOR(10) 
REAL DIFF(10) 
U1 = RMEDIN(X,NWATCH) 
UHAT = Ul 
ALSTOP=0 

C... COMPUTE THE ABSOLUTE VALUE BETWEEN THE DIFFERENCES IN X(I) AND MEDIA 


60 DO 75 I = 1,NWATCH 
IF ( ((X(1)-UHAT)/(6*S))**2 .LT. 1.) THEN 
W=(1-((X(1)-UHAT)/(6*S) ) **2)**2 


ELSE 
W=0 
END IF 
SUM = SUM + W*X(I) 
SUM2 = SUM2 + W 
75 CONTINUE 
TF ( SUM2_LT. .0000001) THEN 
ELSE 
U2=SUM/SUM2 
END IF 
IF(ALSTOP .GT._ 1000)THEN 
TELAG=IFLAG+1 
UHAT=U1 
GO TO 100 
END IF 


LOS 





1F(ABS(U2-UHAT) .GT. .0001) THEN 


AT = U2 
STOP=ALSTOP+1 
TO 10 


Orc 


H 
L 
0 
END IF 
100 BIWAIT=COMPAR(AMEAN, IBIRD, ITRUE,UHAT , BIWAIT ) 
RETURN 
END 
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CEEDS 1000 
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* ITERATIVE REWEIGHTING TO FIND THETA * 


DEFINITION OF LOCAL VARIABLES 


C 
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ese nanos I TRUE) 


Oo 


fe 


)/(SIGMP(1)*DFP))) 


xe 


RMEDIN(Y , NUMB ) 
.NOW,COMPUTE THE VALUE FOR THE WEIGHTS AND THE NEXT THETA 


ETA=TOP/BOTTOM 


iD 


DO 75 I =1,NUMB 


CONTINUE 


. COMPUTE MEDIAN OF OBSERVATIONS 
THETA 
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C2 
Ces 


)"€0 TO 110 


THETA .LT. UP .AND. THETA .GT. LOW 
THETA 


H 
F 
Q 
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T 
I 
H 
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81RD, ITRUE, THETA COUNT) 
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HOLD+. 0001 


LOW=HOLD-. 0001 
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CORO 10 
END IF 


100 COUNT=COMPAR(AMEAN,, IBIRD, ITRUE , UHAT , COUNT) 


RETURN 
END 
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* FUNCTION TO FIND THE MEDIAN * 
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. FINDS THE MEDIAN OF A SET OF NUMBERS 


DEFINITIONS OF LOCAL VARIABLES 


rn 
TUT UTM 
a 


wae 
m 
> 
mo) 
~— bH 


DO 5 
2 CONTENU ie 


_. COMPUTE MEDIAN OF OBSERVATIONS 


DO 25 NM = 1, NUMB 
SMALL=100000. 0 
DO 30 M=1,NUMB 

ie ‘a ee yah) THEN 


END Ik 
30 CONTINUE 
SaaS 
23 CONTINUE 


IF (NUMB/2. .NE. NWYB/2) THEN 
RMEDIN = SMAL TK CNUMB/2. +.5)) 


RMEDIN =(SMAL(NUMB/2)+SMAL((NUMB/2)+1) )/2. 
RETURN 
END 
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* FUNCTION TO COMPARE ESTIMATED MEAN TO TRUE MEAN 
Ke He Fe He He He oe te ke ke te ee KO ke ee 


* 
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... FINDS MEAN WHICH THE ESTIMATE IS CLOSET TO. 
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C DEFINITIONS OF LOCAL VARIABLES 
C COUNT - CURRENT COUNT OF CORRECT ANSWERS 
é ITRUE - TRUE VALUE OF CURRENT BIRDS MEAN 
C MEAN.  - ARRAY CONTAINING ASSUMED MEANS 
¢ NBIRD - INPUTED NUMBER OF BIRDS 
€ UHAT  - CURRENT ESTIMATE OF THE MEAN 
FUNCTION COMPAR( MEAN NBIRDS , ITRUE , UHAT , COUNT ) 
REAL MEAN( 10) ,UHAT, COUNT 
COMPAR=COUNT 
IF ( UHAT .LT. MEAN(1) .AND. MEAN(1) . EQ. ITRUE)THEN 
COMPAR=COUNT+1 
RETURN 
ELSE IF(UHAT. GT. MEAN(NBIRDS). AND. MEAN(NBIRDS). EQ. ITRUE)THEN 
COMPAR= COUNT + 1 
ail RETURN 
IF (UHAT. GT ITRUE-.5 . AND. UHAT .LT. ITRUE+.5) COMPAR=COUNT+1 
RETURN 
bs END 
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APPENDIX C 
SIMULATION CODE FOR THE TWO PARAMETER CASE 
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